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1. Introduction 

In 1918 Schottky discovered that the fluctuations in vacuum diodes can be 
related to the discrete nature of the charge carriers [1] . His observation was 
that the power spectrum of the current fluctuations gave direct access to 
the charge e of the discrete carriers responsible for the current. From his 
theoretical considerations he found a relation between the noise power of 
the current fluctuations Si and the average current /, 

S! = 2el, (1) 

a result nowadays known as the Schottky formula. Its consequence is, that 
the current noise provides information on the transport process, which is 
not accessible through conductance measurements only. 

Studies of the noise properties of tunnel-junctions renewed the interest 
in noise later [2]. Since about ten years the investigation of transport in 
quantum coherent structures has boosted the interest in the theory of cur- 
rent noise in mesoscopic structures [3, 4]. Correlations in the transport of 
fermions have lead to a number of interesting predictions. For example, the 
noise of a single channel quantum contact of transparency T at zero tem- 
perature has the form Si = 2el(l — T) [5, 6]. The noise is thus suppressed 
below the Schottky value, Eq. (1). The suppression is a direct consequence 
of the Pauli principle, and is therefore specific to electrons. Particles with 
bosonic statistics or classical particles show a different behaviour, e. g. for 
Bosons the noise is enhanced by a factor 1 + T. A convenient measure of the 
deviation from the Schottky result is the so-called Fano factor F = Si/2el. 
For a number of generic conductors, it turns out that the suppression of the 
Fano factor is universal, i. e. it does not depend on details of the conduc- 
tor like geometry or impurity concentration. A diffusive metal with purely 
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elastic scattering leads to Fdiff = 1/3 [7, 8], which is independent on the 
concrete shape of the conductor [9]. In a symmetric double tunnel junction, 
on the other hand, F^bitun = 1/2 [10]. A chaotic cavity (a small region with 
classical chaotic dynamics, connected to two leads by open quantum point 
contacts) shows a suppression of F cav i ty = 1/4 [11]. Thus, we conclude that 
from a noise measurement two kinds of information can be obtained. First, 
we can get information on the statistics of the carriers, e. g. if they are 
fermions. Second, provided the statistics is known, the comparison of the 
magnitude of the noise power with the average current gives information on 
the internal structure of the system. However, the picture described here is 
a little bit too simplified. In many real experiments the structure is much 
less well defined, the temperature is finite, or other complications make the 
interpretation of experimental data less trivial. 

Having motivated the interest in the noise, we may ask, what we can 
learn further from current fluctuations. Higher correlators of the current 
will provide even more information on the transport process. However, the- 
oretical calculations of higher correlators become increasingly cumbersome 
and one should find a different concept to obtain this information. This 
step was performed by a transfer of concepts from the field of quantum 
optics. Here it is possible to count experimentally the number of photons 
occupying a certain quantum state. This number is subject to quantum and 
thermal fluctuations and requires a statistical description: the full counting 
statistics (FCS). Lesovik and Levitov adopted this terminology to meso- 
scopic electron transport [12, 13], in which the electrons passing a certain 
conductor are counted. Since then the FCS has been studied in the field of 
mesoscopic electron transport. It the aim of this article to review some of 
the progress, which has been recently made. 

The general problem of counting statistics has been considered before 
on a heuristic level. If one makes the ad hoc assumption that individual 
transfers of charges are uncorrelated and unidirectional, simple calculation 
show that the probability distribution of the number of transfered charges 
is Poissonian 

P^N)=e-*™*&f. (2) 

Here to is the time period during which the charges are counted and N(to) = 
Ito/e is the average number of transfered charges. Schottky's result (1) for 
the current noise power can be easily derived from Eq. (2). This kind of 
counting statistics is usually found in tunnel junctions, where the charge 
transfers are rare events, or at high temperatures, when the mean occupa- 
tion of the states is small and the statistics of the particles plays no role. 
However, in a degenerate electron gas one encounters a completely different 
situation: all states are filled due to Fermi correlations. If we now consider a 
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quantum transport channel of transmission probability T and applied bias 
voltage V, the rigidity of the Fermi sea leads to a fixed number of electrons, 
which are sent into each quantum channel. A charge is transfered to the 
other side with probability T and the statistics is therefore binomial 

p t0 (N) = ^ M ^yN {1 _ T) M(t 0) -N_ (3) 

This is the result of the quantum calculation of Levitov and coworkers [12]. 
Note, that arguments based on the FCS have been used already to inter- 
pret current noise calculations or measurements. However, in general these 
interpretations are not unique, and a calculation of the FCS is required to 
interpret the results for the noise unambiguously. Thus, obtaining the FCS 
is a theoretical task, which leads to a better understanding of quantum 
transport processes. 

The structure of this article is as follows. In the next section we intro- 
duce the theoretical basis, necessary to obtain the FCS. Our approach is 
based on an extension of the well-known Keldysh-Green's function tech- 
nique [14] and described in Section 2.2. A convenient simplification is the 
circuit theory of mesoscopic transport, developed by Nazarov [15], which 
allows to obtain the FCS for a large variety of multi-terminal structures 
with minimal calculational overhead. In the two following sections, we will 
discuss several concrete examples. First, we show that for phase coherent 
two-terminal conductors the counting statistics can be obtained in quite a 
general form [16]. This is illustrated for normal and superconducting con- 
strictions as well as two-barrier structures. A somewhat special case is the 
diffusive conductor coupled to a superconductor (i. e. in the presence of 
proximity effect), where intrinsic decoherence between electrons and holes 
influences the transport properties [17]. Second, we turn to more complex 
mesoscopic structures with more than two terminals. An analytic solution 
is obtained for the case of an arbitrary number of terminal connected by 
tunnel junctions to one central node [18]. As specific example we calculate 
the counting statistics of a beam splitter, in which a normal current or a 
supercurrent is divided into two normal currents. Afterwards we present 
some conclusions. In the appendices we give some theoretical background 
information to the methods used in this article. 

In this article we concentrate on counting statistics in two- or multi- 
terminal devices with either normal- or superconducting contacts. More 
aspects of FCS are in other parts of this book. We would like to mention 
works related to FCS not covered here. Normal-superconductor transport at 
finite energies and magnetic fields was recently addressed in [19, 20]. Time- 
dependent transport phenomena in normal contacts have been studied in 
[13]. Fluctuations of the current in adiabatic quantum pumps have attracted 
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some attention (see e.g. [21]). A connection to photon counting or photon 
transport can be found in Refs. [22, 23]. Counting statistics has also been 
addressed in the context of the readout process of a qubit [24] or to study 
spin coherence effects [25]. Results for the FCS of entangled electron pairs 
[26] and of resonant Cooper pair tunneling [27] have also been published. 



2. Counting Statistics and Green's function 

2.1. BASICS OF CURRENT STATISTICS 

We introduce some basic formulas, relevant for the theory of FCS. The 
quantity we are after is the probability Pt (N), that N charges are trans- 
fered in the time interval to- Equivalently, we can find the cumulant gener- 
ating function (CGF) S(x), defined by 

exp(-S( x )) = E P to( N ) exp(iiV x ) . (4) 

N 

To keep notations simple in this section, we will limit the discussion to the 
two-terminal case, in which only the number N of transfered charges in one 
terminal matters. In the other terminal it is given by — N, since the total 
number of charges is conserved. However, most relations are straightfor- 
wardly generalized to many terminals. Note that normalization of the FCS 
requires that 5(0) = 0. Also, we will suppress the explicit dependence of 
S(x) on the measuring time to- In the static case considered mostly in this 
article, we have S(x) ~ to- 

From the full counting statistics one obtains the cumulants 



C 1 = N^J2 NP to( N ) > C 2 = (N-N)\ (5) 

N 



=—2 



C 3 = (N - N) 3 , C 4 = (N - N) 4 - 3{N - N) 2 , (6) 

and so on. The meanings of the various cumulants are depicted in Fig. la. 
Alternatively, the cumulants can be obtained from the CGF 



(7) 



The relation to the average current and the noise power of current fluctu- 
ations is obtained as follows. Almost trivially one has 

C 1 = (N) = — [ t0 dt(I(t)) = -^I, (8) 
-e Jo e 
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where we have denoted the charge of the electrons with — e. The relation 
between the second cumulant and the current noise power 



Si 



/oo 
dr({SI(r),SI(0)}} (9) 



is less obvious. We write for the second cumulant 

C 2 = (A -AT) 2 = ± jT° jT° dtdt'({SI(t), (57(f)}) , (10) 

where SI(t) = I(t) — (I) is the current fluctuation operator and (...) denote 
the quantum statistical average. We transform the time coordinates to av- 
erage t = (t + t')/2 and relative time r = t — t'. Assuming the observation 
to is much larger than the correlation time of the currents, the correlator 
in Eq. (10) does not depend on T, and we find the desired relation between 
current noise power and the second cumulant 



Si 2e 2 d*S( X ) 



to dx 2 

For higher correlators similar formulas can be derived 



(11) 



2.2. EXTENDED KELDYSH GREEN'S-FUNCTION TECHNIQUE 

The task to measure the number of charges transfered in a quantum trans- 
port process has, in general, to be formulated as a quantum measurement 
problem. A thorough derivation goes beyond the scope of this article and 
we refer to Ref. [28]. The quantum- mechanical form of the cumulant gen- 
erating function is given by [14, 17, 28] 

e -S(x) = (Te^ Io° dt mf e -i£ C dm)) . (12) 

Here T(T) denotes (anti-)time ordering and I(t) the operator of the current 
through a certain cross section. As preliminary justification of Eq. (12) 
we note, that it is easily shown that expansion of Eq. (12) in % yields 
the various current-correlators. The expectation value in Eq. (12) can be 
implemented on the Keldysh contour, see Appendix A.l, which makes it 
possible to use diagrammatic methods [29]. Equation (12) has a form similar 
to the thermodynamic potential in an external field. From the linked cluster 
theorem (see Appendix A. 2) it follows that the CGF is the sum of all 
connected diagrams. 

To connect the CGF to accessible field-theoretical quantities we con- 
sider the nonlinear response of our electronic circuit to the time-dependent 
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Figure 1. Left panel: An example of a probability distribution illustrating the meaning 
of the different cumulants. The average is given by Ci, the width by C2, the 'skewness' 
by C3 and the 'sharpness' is related to the C4. Right panel: schematics of an ideal charge 
counter in a terminal. The number of charges passing the cross section C are counted. 



perturbation 

H c (t) = |-I c (i) = T |- J tfW(x, i)j c (x)*(x, t) , (13) 

where the +(— ) sign is taken on the lower (upper) part of the Keldysh con- 
tour. The operator j c ( x ) is the operator of the current through a cross 
section c, depicted in Fig. 1. We allow here for multi-component field op- 
erators, such as spin or Nambu for example. Matrices in this subspace are 
denoted with a Since we are aiming at the total charge counting statis- 
tics, we will assume that x is nonzero and constant in a finite time interval 
[0,to]. 

The unperturbed system evolves according to a Hamiltonian Ho = 
J d?x^(x)ho(x)ty(x), where ho is the usual single-particle Hamiltonian of 
the system. The equation of motion for the matrix Green's function subject 
to H = Hq + H c (t) reads (in the Keldysh matrix representation) 



d * y » 



Gtx^x'^xH^-O^x-x), (14) 



Here T3 denotes the third Pauli matrix in the Keldysh space. The relation of 
the Green's function (14) to the CGF (12) is obtained from a diagrammatic 
expansion in x- One finds the simple relation (see App. A. 2) 

^ = ^I(X). (15) 
ox e 

The counting current I(x) is obtained from the x-dependent Green's func- 
tion via 

I( X ) = I d 3 x Tr [r 3 j c (x)G(x, t; x', t; X )l , . (16) 

Since we are assuming a static situation, the r.h.s. of Eq. (16) is time- 
independent. The relations (14)-(16) offer a very general way to obtain the 
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full counting statistics of any system. It is, however, difficult to find the 
Green's function in the general case. 

For a mesoscopic transport problem there is a particularly simple way to 
access the full counting statistics, based on the separation into terminals (or 
reservoirs) and an active part, the first providing boundary conditions, and 
the second being responsible for the resistance. Let us consider the equa- 
tion of motion for the Green's function inside a terminal for the following 
parameterization of the current operator in Eq. (13) 

if* 

j c (x) = (VF(x)) lim — (V x - V*) <7 3 . (17) 

F(x) is chosen such that it changes from to 1 across a cross section C, 
which intersects the terminal, but is of arbitrary shape, see Fig. 1. Here we 
have added a matrix 03 to the current operator, which accounts for possible 
multicomponent field operators like in the case of superconductivity. The 
change from to 1 should occur on a length scale A, for which we assume 
Xf <C A <C limp,^o (Fermi wave length Xp, impurity mean free path li mp , 
and coherence length £0 = vp/2A). Under this assumption we can reduce 
Eq. (14) to its quasiclassical version (see Ref. [30] and App. B.l). This is 
usually a very good approximation, since all currents in a real experiment 
are measured in normal Fermi- liquid leads. 

The Eilenberger equation in the vicinity of the cross section reads 



v F V<7(x,v F ,M',x) = 



-i-(VF(x))v F fx , <7(x,v F ,M',x) 



(18) 



Here fx = T3CT3 is the matrix of the current operator. Other terms can be 
neglected due to the assumptions we have made for A. The counting field 
can then be eliminated by the gauge-like transformation 

S(x, v P , t, t', X ) = e-^ F ^l 2 g(^ v F , t, t', yxF(x)f K /2 

We assume now that the terminal is a diffusive metal of negligible resistance. 
Then the Green's functions are constant in space (except in the vicinity 
of the cross section C) and isotropic in momentum space. Applying the 
diffusive approximation [31] in the terminal leads to a transformed terminal 
Green's function 

G( x ) = e- ixfK/2 G(0)e ixfK/2 , (20) 

on the right of the cross section C (where F(pc = 1) with respect to the case 
without counting field. Consequently the counting field is entirely incorpo- 
rated into a modified boundary condition imposed by the terminal onto the 
mesoscopic system. Note, that it follows from (17) and (18), that the count- 
ing field for a particular terminal vanishes from the equations of motion in 
the mesoscopic system and the other terminals. 
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The generalization of this method to the counting statistics for multi- 
terminal structures was performed in [32]. The surprisingly simple result 
is, that one has to add a separate counting field for each terminal, in which 
charges are counted. In Sec. 4.1 we demonstrate this for an example, in 
which an analytic solution can be found. This concludes the derivation of 
our theoretical method to obtain the FCS. 

What are the achievements of this method? We should emphasize that 
it does not simplify the solution of a specific transport problem, i.e. we still 
have to know the solution corresponding to the Hamiltonian Hq. If this 
solution is not known, the counting field makes this situation not easier. 
Rather, the method paves a very general way to obtain the FCS, if a method 
to find the average currents, i.e. for x = 0, is known. In the next section we 
will introduce such a method, the circuit theory of mesoscopic transport. 
Initially it was invented to calculate average currents, however the method 
to obtain the FCS introduced in this section is straightforwardly included. 

What is the price to pay? Loosely speaking, the method to obtain the 
average currents has to be sufficiently general. Usually the absence of a 
field, which has different signs on the upper and lower part of the Keldysh 
contour, allows some simplification. For example, in the Keldysh-matrix 
representation all Green's functions can be brought into a tri-diagonal form, 
which is obviously simpler to handle than the full matrix. The method above 
does not allow this simplification anymore. Or, in other words, the counting 
rotation (20) destroys the triangular form. Thus, the price we have to pay 
for an easy determination of the FCS is that we need a method, which 
respects the full Keldysh-matrix structure in all steps. The circuit theory, 
which we describe in the next section fulfills this requirement. 

2.3. CIRCUIT THEORY 

A concise formulation of mesoscopic transport is the so-called circuit theory 
[15, 33]. Its main idea, borrowed from Kirchhoff's classical circuit theory, is 
to represent a mesoscopic device by discrete elements. These approximate 
the layout of an experimental device with arbitrary accuracy, provided one 
chooses enough elements. In practice one has to find the balance between 
a small grid size and the computational effort. 

We briefly repeat the essentials of the circuit theory. Topologically, one 
distinguishes three elements: terminals, nodes and connectors. Terminals 
are the connections to the external measurement circuit and provide bound- 
ary conditions, specifying externally applied voltages, currents or phase dif- 
ferences. Besides, they also determine the type of the terminal, i.e. if it is a 
normal metal or a superconductor. The actual circuit, which is to be stud- 
ied consists of nodes and connectors, the first determining the approximate 



9 



layout, and the second describing the connections between different nodes. 

The central element of the circuit theory is the arbitrary connector, 
characterized by a set of transmission coefficients {T n }. Its transport prop- 
erties are described by a matrix current found in [33] 



e 



2 _ 2T, 



n , G2 



112 = — E / - - — v- ( 21 ) 

* n 4 + T n ({Gi,G 2 }-2) 

Here G , 1 ( 2 ) denote the matrix Green's functions on the left and the right of 
the contact. We should emphasize that the matrix form of (21) is crucial 
to obtain the FCS, since it is valid for any matrix structure of the Green's 
functions. The electrical current is obtained from the matrix current by 



'12 



±- J dETrf K I 12 . (22) 



A special case is a diffusive wire of length L in the presence of proximity 
effect. If L is longer than other characteristic lengths like £0, decoherence 
between electrons and holes becomes important and the transmission eigen- 
value ensemble is not known. Instead, one solves the diffusion-like equation 
[31] 

VD(x)G(x)VG(x) = \-iEf 3 , G(x)l . (23) 

The matrix current is now given by 

/(x) = -(j(x)G(x)VG(x) . (24) 

In these equations, -D(x) the diffusion coefficient and cr(x) = 2e 2 NqD(x.) is 
the conductivity. In general this equation can only be solved numerically, 
but in some special cases (e.g. for E = 0) an analytic solution is possible. 

If the circuit consists of more than one connector, the transport proper- 
ties can be found from the circuit theory by means of the following circuit 
rules. We take the Green's functions of the terminals as given and intro- 
duce for each internal node an (unknown) Green's function. The two rules 
determining the transport properties of the circuit completely are 

1. G 2 j = 1 for the Green's functions of all internal nodes j. 

2. The total matrix current in a node is conserved: J2i hj = 0, where the 
sum goes over all nodes or terminals connected to node j and each 
matrix current is given by (21) or (24), depending on the type of the 
connector. 

An important feature of the circuit theory in the form presented above is 
that it accounts for any matrix structure (i.e. Keldysh, Nambu, Spin, etc.). 
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Figure 2. Two terminal contacts. The leftmost picture shows a general two-terminal 
contact. A matrix current I is flowing between terminals, characterized by Green's func- 
tions Gi,2- Several connectors (a-e) are described in the text and depicted as: a) tunnel 
junction, b) arbitrary connector with transmission eigenvalues {T n }, c) two tunnel junc- 
tions in series with conductances gjv(ss), respectively, d) a chaotic cavity connected by 
N c h channels to the leads, and e) a diffusive wire, characterized by a resistance Rn- 



Thus, we can straightforwardly obtain the FCS along the lines of Sec. 2.2. If 
the charges in a terminal are counted, we have to apply a counting rotation 
(20) to the terminal Green's function. The counting-rotation matrix has the 
form fx = 03T3, where o"j(fj) denote Pauli-matrices in Nambu(Keldysh)- 
space. Then we can employ the circuit rules to find the x-dependent Green's 
function and finally obtain the total CGF by integrating all currents into 
the terminals over their respective counting fields (see [32] and [18] for more 
details). 



3. Two- Terminal Contacts 

In this chapter we demonstrate several examples of the FCS of contacts 
between two terminals, see Fig. 2. One can easily derive a number of rather 
general results, such as Poisson statistics in the case of tunnel junction, 
or binomial statistics for single channel contacts of transparency T. All 
these results can be found from a general CGF [16], which depends on the 
ensemble of transmission eigenvalues {T n }. To illustrate specific examples, 
we will compare the cases of transport between two normal terminals or 
between one superconducting and one normal terminal. In the end, the 
statistics of an equilibrium supercurrent is discussed. 



3.1. TUNNEL JUNCTION 



The counting statistics of a tunnel junction contact can be obtained from a 
direct expansion of the matrix current (21), if T n <C 1 for all n. It coincides 
with the result obtained by means of the tunneling Hamiltonian [34] . The 
matrix current takes the form 



(25) 



Here the matrix current depends only on the tunneling conductance = 
(e 2 /-7r) J2 n Tn of the contact, and we have arbitrarily chosen to count the 
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charges in terminal 1. Now, the counting current is 



7 (x) = -gf / dETr{f K [G 1 (x),G 2 \) . 



(26) 



We use the pseudo-unitarity f\ = 1 to express the counting rotation as 
exp(ixTK-/2) = f+exp(ix/2) + f_ exp(— z%/2), where f± = (1 ± tk)/2. 
Then the current has the form 



J(x) = 



to 
tpGr 
4e 2 



21 e 



/ 



dETr 



f + Gif-Gj 



(27) 
(28) 



The CGF follows from integrating (27) with respect to x an d we obtain 

S(x) = -N 12 (<? x - 1) - JV 2 i(e _i * - 1) . (29) 

It is easy to see that the even and odd cumulants obey 

C 2 „+i = N 12 - N 21 , C 2n = N 12 + N 21 . (30) 

If only tunneling processes in one direction occur (say from 1 to 2), A^i = 
and the average N = Ni 2 . The statistics is Poissonian 



Pto(N) 



-N 



N 



N 



Nl 



(31) 



In particular, for the current noise we find the Schottky result Si = 2el. 

We conclude that the charge counting statistics of a tunnel junction 
(or more precisely, if the transfer events are rare) is of a generalized Pois- 
son type [34]. If only tunneling events in one directions are possible, the 
statistics is Poissonian. 



3.2. QUANTUM CONTACT - GENERAL CONNECTOR 

We consider now a quantum contact characterized by a set of transmission 
eigenvalues {T n }. It turns out that the CGF can be obtained in a quite 
general form, valid for arbitrary junctions between superconductors and/or 
normal metals. 

The matrix current through a quantum contact is described by Eq. (21) 
and the CGF can then be found from the relation dS(x)/°X = ( — ^o/4e 2 ) 
/ dE Tr[f K I(x)]. Using that [A,{A,G 2 }} = for all matrices with A 2 = 1 
the following identity holds 

Tr^-{G 1 ( X ),G 2 ] n = l -Trf K [(5i(x),G< 2 ] {G^x), G^* . (32) 
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We can therefore integrate (22) with respect to \ and obtain [16] 



(33) 



This is a very important result. It shows that the counting statistics of 
a large class of constrictions can be cast in a common form, independent 
of the contact types. Another important property of Eq. (33) is, that the 
CGF's of all constrictions are linear statistics of the transmission eigenvalue 
distribution, and can therefore be averaged over by standard means (see 
e.g. [35]). Examples are given in the following sections. Note also, that an 
expansion of Eq. (33) for T„<Cl yields the result (29) of the tunnel case. 

We will now discuss several illustrative examples. Consider first two 
normal reservoirs with occupation factors f\( 2 ) = [ ex P((^ — Mi(2))/^e) + l] _1 
(T e is the temperature). We obtain the result [12, 13] 



S M = ~ £/^hi [l + B ln (E) (e* - l) + B. ln (E) (e-* - l) 



.(34) 



Here we introduced the probabilities B\ n = T n f\(E) (1 — /2(E)) for a tun- 
neling event from 1 to 2 and B_\ n for the reverse process. The terms with 
counting factors e ±tx — 1 obviously correspond to charge transfers from 1 
to 2 (2 to 1). At zero temperature and — Hi = eV > the integration 
can easily be evaluated and we obtain 



Six) 



et Q \V\ 



TT 



[l + T n (e*-l 



(35) 



The corresponding statistics for a single channel with transparency T is 
binomial 



PtoiN) 



( ^\T N (I — T) M ~ ' 



(36) 



Here we have introduced the number of attempts M = etoV/w, which 
is the maximal number of electrons that can be sent through one (spin- 
degenerate) channel in a time interval to due to the exclusion principle for 
Fermions. 

The FCS of an superconductor-normal metal-contact also follows from 
Eq. (33) (for a definition of the various Green's functions see App. B.4). 
Evaluating the trace in Eq. (33) the CGF can be expressed as [36] 



six) 



dE In 



1+ £ A nq (E) (e^ - l) 



(37) 
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The coefficients A nq (E) are related to a charge transfer of q x e. For example 
a term exp(2i%) — 1 corresponds just to an Andreev reflection process, in 
which two charges are transfered simultaneously [37]. Explicit expressions 
for the various coefficients are given in Appendix C.l. For the BCS case, 
they reproduce the results of Ref. [36]. In the fully gapped single-channel 
case at energies ksT e <C eV <C A only terms corresponding to Andreev 
reflection (A2) are nonzero and the CGF becomes 



six) 



et \V\ 



In 



TT 



1 + 



Ra (t 



i2 X 



1 



(38) 



where Ra = T 2 /(2 - T) 2 is the probability of Andreev reflection. The CGF 
is now 7r-periodic, which means that only even numbers of charges are trans- 
fered, a consequence of Andreev reflection. The corresponding statistics is 
binomial 



Pt m) = 



\R%{l-R A ) M - N ,P to (2N+l) = 0. 



(39) 



The number of attempts M is, however, the same as in the normal state. 

It is interesting to see how the CGF for normal transport, i. e., Eq. (34), 
emerges from Eq. (37). Putting fji jA = and gR jA = ±1 the coefficients in 
Eq. (37) can be written as 



A ±2 



±1 



B+ 1 + B 

B ±i B ±\ ; 



± 



Bf = TBi(±E) ,B± X = TB_i(±£?) , 



The argument of the 'In' in Eq. (37) factorizes in positive and negative 
energy contributions 



In 



(40) 



Integrating over energy both terms give the same contribution, and the 
CGF results in Eq. (34). This shows explicitly how the positive and negative 
energy quasiparticles are correlated in the Andreev reflection process. 



3.3. DOUBLE TUNNEL JUNCTION 

We now consider a diffusive island (or a chaotic cavity) connected to two 
terminals by tunnel junctions with respective conductance g± and 52 [18]. 
We assume for the conductance of the island gisiand S> 51,2 S> e 2 /h, so 
we can neglect charging effects. This provides a simple application of the 
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circuit theory. The layout is shown in Fig. 2c. The central node is described 
by an unknown Green's function G c . We have two matrix currents entering 
the node, which obey a conservation law: 



= h + h = j 



giG l +g 2 G 2 ,G c \ . (41) 
Using the normalization condition G 2 C = 1 the solution is 



G c = , 9iG 1+ g 2 G 2 (42) 



hi + ^ +01S2 {6i,<3 2 } 
We can integrate the current /(%) ~ Trfxli and obtain the CGF 

Six) = -^J dETr s jgl + gl+g 1 g 2 {G 1: G 2 } . (43) 

Again, this result is valid for all types of contacts between normal metals 
and superconductors. 

We first evaluate the trace for two normal leads and find 

S( X ) = - t fJdEx (44) 

yj(9i +g2? +4<?i 52 (h(l - f 2 )(e*x - 1) + f 2 (l - h)(e^x _ 1)) . 

We observe that the CGF contains again counting factors corresponding to 
charge transfer from 1 to 2 and vice versa. In contrast to the Poissonian case 
for a tunnel junction, Eq. (29), the charge transfers are not independent, but 
correlated by the square-root function. At zero temperature and fi± — fi 2 = 
eV > we find the result, 

s(x) = -^y/ioi + g2? + ±gi 9 2(e*x - 1) . (45) 

There are two relatively simple limits. If the two conductances are very 
different (e.g. g\ <C g 2 ), we return to Poissonian statistics: 

S(x) = -^(e*-l). (46) 
On the other hand, in the symmetric case g\ = g 2 = g we find [38] 

5(x) = -^(e^ 2 -l), (47) 
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and the cumulants are 

C- * , N = ^. (48 ) 

The suppression factor 1/2, which occurs already in the Fano factor, carries 
forward to all cumulants. Note, that the same kind of statistics (44) follows 
also from a Master equation [38, 39]. 

The CGF (43) for the transport between a normal metal and a super- 
conductor (at zero temperature, for simplicity) reads [18] 



S(x) = J ^al + al + \]{al + alf + ^(e** - l) . (49) 

Thus, the influence of the superconductor is two-fold: charges are transfered 
in units of 2e (indicated by the 7r-periodicity) and another square root is 
involved in the CGF, resulting from the higher order correlations. In the 
limit that both conductances are very different (e.g. g\ <C 52) we obtain 
again Poissonian statistics 



*(X) = -^1(^-1) • (50) 



This corresponds to uncorrelated transfers of pairs of charges. Consequently 
we obtain for the cumulants 

C n = 2 n ~ 1 N , N=^^, (51) 

and the effective charge 2e can indeed be found from the Fano factor. The 
transport properties at finite energies and magnetic fields of this structure 
have recently been addressed in [19] and [20]. 

3.4. SYMMETRIC CHAOTIC CAVITY 

Another interesting system is the chaotic cavity, i. e. a small island coupled 
to terminals by perfectly transmitting contacts (with N c f- L 3> 1 channels). 
This system is described by the circuit depicted in Fig. 2d. The matrix 
current between terminal 1 and the cavity is I\ = N^e /tt)[Gi, G c ]/(2 + 
{Gi, G c }). Similar as in the previous chapter, the current conservation reads 
now 

G4,GU G2,GU 

= L ( z I 1 + L , z , (52) 

2 + {Gi,G c } 2 + {G 2 ,G c } 



16 



which is solved by the solution (42) for g\ = g2- The integration of (22) 
leads to 



Six) = ~ N ch^ J dETr In 2 + ^2 + {(5i(x) , G 2 } 



(53) 



The interpretation of this result is straightforward. As we have seen in 
Sec. 3.2, the In appeared already in the FCS of a quantum point contact 
(leading to binomial statistics). The square-root we encountered already in 
the previous section and we attribute it to inter-mode mixing on the central 
node ('cavity noise'). 

For normal leads at zero temperature and applied bias voltage V we 
obtain (with the number of attempts M = N c f l toeV/'K) 



S(x) = -2Mln 



1 + e 



On the other hand, in the case of Andreev reflection we find 



Six) = -2Mln 



2 + e *[xmod7r] + + e i[xmod7r] 



(54) 



(55) 



where the 7r-periodicity reflects the fact that charges are transfered in pairs. 



3.5. DIFFUSIVE CONNECTOR 

A metallic strip of length L with purely elastic scattering, characterized 
by an elastic mean free path I, is called a diffusive connector if Z <C L. 
Its transport properties are governed by the diffusion-like Usadel equation, 
Eq. (23). In the case of proximity transport the r.h.s. of Eq. (23) accounts 
for decoherence of electrons and hole during their diffusive motion along 
the normal wire. This term has the form of a leakage current, if the l.h.s. is 
considered as a conservation law for the matrix current [33]. Note that the 
electric current is still conserved, it is only loss of coherence which occurs. 
In general the solution of the full equation is rather complicated and can 
only be found numerically for the full parameter range. There is, however, 
one case, in which an analytic solution is possible, namely if the r.h.s. of 
Eq. (23) vanishes. This is either the case for purely normal transport, when 
electrons and holes are transported independently, or for E = 0, which 
means we are restricted to low temperatures and voltages. The scale here 
is set by the Thouless energy Eth (given by %D / L 2 for a wire of uniform 
cross section) . At this scale the famous reentrance effect of the conductance 
occurs [40]. This regime was studied in Ref. [17] and will be discussed in 
connection with experimental results for the current noise in the article 
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by Reulet et al. in this book. Numerical results for equilibrium counting 
statistics in the full parameter range are discussed below. 

We now concentrate on the analytic solution in a quasi-one-dimensional 
geometry, i. e. we assume a wire of uniform cross section connects two 
reservoirs located at x = and x = L. It is characterized by a conductivity 
a(x), which in general could depend on x, e. g. due to an inhomogeneous 
concentration of scattering centers. The diffusion equation is then indeed a 
conservation law for the matrix current density 

JLj = , j = -a(x)G(x)-^-G(x). (56) 

This equation has to be solved with the boundary condition that G(0) = G\ 
and G(L) = Gi- It follows from Eq. (56) and the normalization condition 
G 2 (x) = 1 that G(x) obeys the equation 

d ~ 1 - 

j = const. , -k-G{x) = -—G(x)j . (57) 

ox a(x) 

This is an homogeneous first-order differential equation, which is easily 
solved. Using the boundary conditions we find the solution 



G(L) = G(0)e / /^ , g d = " , (58) 

Jo dx/a(x) 

where gd is the conductance of the wire and A its cross section. The current 
is thus given by [14] 

I(x) = -9dln(G 1 ( X )G 2 ) , (59) 

where we have reinserted the dependence on the counting field %. To find 
the CGF we have to find the integral / dxT^fxIix) with respect to the 
counting field. Expanding the In and using repeatedly the normalization 
condition the results is 



A 



S W = - tj B J dETr\ [ acosh 2 (^{G 1 (x),G 2 } 



(60) 



This is the counting statistics for a general diffusive contact (under the 
restrictions mentioned above). The first thing to note is that from the 
properties of the wire only the conductance enters and this holds for all 
cumulants. In this sense Eq. (60) shows that the entire FCS is universal. 
In our derivation, we have assumed a wire of uniform cross section, but it 
has been shown [9] that this also holds for an arbitrary shape of the wire 
(as long as it can be considered as quasi one-dimensional). We should also 
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mention, that one could have obtained the same result, by averaging the 
CGF (12) over the bimodal distribution of transmission eigenvalues [41, 9]. 

Now we evaluate the trace in Eq. (60) for normal metals at zero tem- 
perature and applied bias voltage eV. We obtain 

S (x) = -^g^acosh 2 (2e* - l) , (61) 

which coincides with the results of Refs. [14, 42]. 

In the case of Andreev transport the easiest way to obtain the CGF is 
as follows. We have already previously noted, that Eq. (61) follows from 
averaging Eq. (35) with the transmission eigenvalue distribution p(T) = 
(2e 2 /g d 7t)/TVl-T for a diffusive metal [9, 41]. Now, the CGF for Andreev 
transport Eq. (38) has the same form as in the case of normal transport, 
provided we replace x with 2\ and the transmission eigenvalue T n with the 
Andreev reflection probability Ra(T u ). A simple calculation shows, that 
the Ra are distributed according to the same distribution as the normal 
transmission eigenvalues (up to a factor of 1/2). Thus, we can immediatly 
read off the CGF for the diffusive SN-wire in the limit of zero temperature 
and eV <C E^h from Eq. (61) and obtain 

S( x ) = -^^acosh 2 (2e 2i * - l) . (62) 
8e ^ ' 

As a consequence the relation between the cumulants in the SN-case, C^ N , 
and the normal case, C^ N , is 

= 2"- 1 C^ N . (63) 

We observe that we can read off the effective charge from the ratio C^ N /C^ N 
= (q e ff/e) n ~ 1 and, indeed, find q e ff = 2e. We should emphasize, how- 
ever, that this is a special property of the diffusive connector. Our prove of 
Eq. (63) is valid as long as eV <C Ej>h, and it is not clear, wether Eq. (63) 
is true also for eV S> E^h 

The counting statistics in equilibrium for arbitrary temperature T e was 
studied in Ref. [17]. By a numerical solution of Eq. (23), it is possible 
to evaluate the integral over x i n the inversion of Eq. (4) in the saddle 
point approximation, i.e. we take % as complex and expand the expo- 
nent around the complex saddle point x = ^ x o- The integral yields then 
P(N) « exp(— S(ixo) — xqN), which we plot implicitly as a function of 
N(xq) = dS(ixo)/dxo. Results of this calculation are displayed in Fig. 3. 
The charge number N is normalized by A^h er m = 9dtokBT e /e 2 . Note that 
the conductance (and due to the fluctuation-dissipation theorem the noise) 
is the same in all cases. The solid line shows the distribution in the normal 
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Figure 3. Equilibrium distribution of the current fluctuations in a diffusive SN-contact. 
We observe a) a strong deviation of all distributions from a Gaussian b) enhanced fluc- 
tuations in the superconducting case, and c) differences between proximity effect and 
coherent transport. Note, that the second cumulant (i.e. the thermal noise) is the same 
for all displayed curves. 

state, which does not depend on the Thouless energy. In our units, this curve 
is consequently independent on temperature. In the superconducting state 
the Thouless energy does matter, and the distribution depends on the ratio 
ETh/ksTe. We observe that large fluctuations of the current in the super- 
conducting case are enhanced in comparison to the normal case, and in both 
cases are enhanced in comparison to Gaussian noise ~ exp(— A r2 /4iV t i iem ). 
The differences between the normal and the superconducting state occur 
in the regime of non-Gaussian fluctuations. 

3.6. SUPERCURRENT 

The CGF of a quantum contact, i. e. Eq. (12), can be used to find the 
counting statistics for a supercurrent between two superconductors at a 
fixed phase difference </>. This has been done in Ref. [16]. The result can be 
represented in a form similar to the CGF of the SN-contact (37) 

Explicit expressions for the coefficients are given in Appendix C.2. To find 
the statistics of the charge transfer, we will treat two separate cases. 

Gapped superconductors. If the two leads are gapped like BCS supercon- 
ductors, the spectral functions are given by Eq. (115). Here we account for 
a finite lifetime S of the Andreev bound states, e. g. due to phonon scatter- 
ing. The supercurrent in a one-channel contact of transparency T\ is solely 
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carried by Andreev bound states with energies ±A(1 — T\ sin 2 4>/2) 1 / 2 = 
±Eb ((f))- The importance of these bound states can be seen from the coef- 
ficient Q(E,(fi) (see Eq. (124)), which may become zero and will thus pro- 
duce singularities in the CGF. 1 The broadening 5 shifts the singularities of 
Q(E, 4>) into the complex plane and allows an expansion of the coefficients 
Af close to the bound state energy. Performing the energy integration the 
CGF results in 



S( X , 4>) = -2t SJ I-X'^r- tanh j , 



(65) 



where I\(<p) = A 2 Ti sui((/>)/2.Eb(<^) is the supercurrent carried by one 
bound state. In deriving (65) we have also assumed that % <S 1. This 
corresponds to a restriction to "long trains" of electrons transfered, and 
the discreteness of the electron transfer plays no role here. Fast switching 
events become less probable at low temperatures and are neglected here. 
In the saddle point approximation and for 7(0) = 1/ cosh(£ , s(0)/2T e ) <C 1 
the FCS can be found. We express the transfered charge in terms of the cur- 
rent normalized to the zero temperature supercurrent: j((f>) = eN/toh(4>)- 
We find for the current distribution in the saddle point approximation 



P(j, 4>) ~ - exp 
7 



) , (66) 



for |j(i^>)| < 1 and zero otherwise. At zero temperature Eq. (66) approaches 
P(j,(p) — ► S(j — 1). Thus the charge transfer is noiseless. At finite tem- 
perature, on the other hand, the distribution (66) confirms the picture of 
switching between Andreev states which carry currents in opposite direc- 
tions, suggested in Ref. [43]. The previous result is valid under the following 
conditions. In the energy integration it was assumed that the bound states 
are well defined. For small transmission the distance of the bound state to 
the gap edge is ~ T±A. Thus, to have well-defined bound states we have to 
require S < T\A. Similarly, for a highly transmissive contact and a small 
phase difference we require (f> ~ eN/toI c > S/A. The statistics beyond these 
limits is similar to what is discussed in the following. 

Tunnel junction/ 'gapless superconductors. Let us now consider the su- 
percurrent statistics between two weak superconductors, where the Green's 
functions can be expanded in A for all energies. One can see that this is 
equivalent to the tunnel result (29). It has the form 

SUA) = ~N + {m x - 1) - iMMe"* " 1) • (67) 

1 It is interesting to note that in the limit 8 — » the CGF has poles for energies 
E%{x) = A 2 (l -Tsin 2 (^rp))- The counting field therefore couples directly to the 
phase sensitivity of the Andreev bound states. 



21 



where 

JV±(0 = |(P B (#±W)) ■ (68) 

This form of the CGF shows that the FCS is expressed in terms of super- 
current I s ((f>) and noise P s (4>) only. Supercurrent and current noise are 

I s (<f>) = -^Re J dETr{a 3 [R 1 (E,4>),R 2 (E)]}t&nh(y^) , (69) 
P s (cf>) = -^Re ( dETr {<J 3 ii(£, cl>)a 3 R 2 (E)} 1 - . (70) 

COsh \2k B T e ) 

Here Gt is the normal-state conductance of the contact. Eq. (70) shows 
that P s vanishes at zero temperature, whereas I s vanishes at T c . Therefore, 
there is some crossover temperature below which P s < \I S \. In this limit 
one of the coefficients N± becomes negative and the interpretation, that 
the CGF (67) corresponds to a generalized Poisson distribution, makes no 
sense anymore. In fact, the CGF leads to 'negative probabilities' and does 
not correspond to any probability distribution, The origin of this failure is 
the broken U(l)-symmetry the superconducting state. Nevertheless the FCS 
can be used to predict the outcome of any charge transfer measurement, as 
is discussed in detail in Refs. [16, 28, 44]. 

4. Multi- Terminal Structures 

Many mesoscopic transport experiments are performed in multi-terminal 
configurations. An example is shown in Fig. 4. Due to the quantum nature 
of the charge carriers, interesting non-local effects can appear, such as sen- 
sitivity of measured voltage differences to changes in the setup outside the 
current path. Obviously, the same is true for current fluctuations and for 
the full counting statistics. These are sensitive to the quantum correlations 
between the charge carriers, which can have a nonlocal character. In terms 
of counting statistics this means, for example, that the joint probability to 
count Ni particles in terminal 1 and N 2 particles in terminal 2 can not be 
factorized into separate probabilities for the two events. 

4.1. GENERAL RESULT FOR MULTI-TUNNEL GEOMETRY 

The generalization of the method, introduced in Section 2 for two terminals, 
to many terminals is straightforward [32]. The counting field \ is replaced 
by a vector %, with dimension equal to the number of terminals. For brevity 
we collect the charges passing each terminal into a vector N. The current in 
each terminal is coupled to the respective component of the counting field 
by an expression like (13). Following the procedure outlined in Sections 2.2 
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Figure 4- a) An example for a multi-terminal structure. K terminals are connected to 
a central node by tunnel junctions b) Three-terminal structure. Two different voltage 
configurations are considered. In the first case, all terminals are normal metals. A bias 
voltage is applied between terminal I and terminal 2, whereas terminal 3 is operated 
as a voltage probe (no mean current). In the second case we consider a beam splitter 
configuration. A supercurrent or normal current in terminal 3 is divided in (or merged 
from) two normal currents from terminals I and 2. Here, we assume that the same 
potential is applied to terminals I and 2. 



and 2.3 the result is, that the Green's function of each terminal acquires its 
own counting field %n- The rules, that determine the transport properties, 
remain essentially unchanged. 

The procedure outlined above is best illustrated by an example. We 
consider a node connected to K terminals via tunnel junctions with con- 
ductances gt- This setup is shown in Fig. 4a. Each terminal is described 
by a Green's function Gk{Xk) (k = 1..K), which is related to the termi- 
nal's usual Green's function by a counting rotation (20). We do not need 
to specify yet, whether the terminals are superconducting or normal. The 
goal is to find (for arbitrary applied voltages and temperatures) the joint 
probability Pt (N\, N2, ...,Nk), that N± particles enter through terminal 1, 
N2 particles through terminal 2, and Nk particles through terminal K. 
Correspondingly we define the cumulant generating function 

e-*(x)= J2 Pt (N)e^. (71) 

N lt N2,-,N K 

The central node is described by a Green's function G c (x), which has to be 
determined from the circuit rules. The matrix currents through terminal k 
is given by Ik = — (gk/2)[Gk(Xk),G c (x)] and current conservation on the 
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node can be written as 



K 1 



k=l 



K 



k=l 



(72) 



This equation (together with the normalization condition G c {x) = 1) is 
solved by 



G C (X) 



Ek=i9kG k (Xk) 



T,k,i=i9k9i {G k (Xk),Gi{xi)} 



(73) 



The CGF is found from integrating the relations dS(x) / d\k = (—ito/e) 
I k (x), where I k (x) = (V 4e ) / dETrf K I k (x). The CGF is then determined 
up to an additive constant, which is fixed by the normalization S(0) = 
and neglected in the following. The result is [18] 



S(X) = ~^jdE T^Y.l l=l 9k9i{Gk{Xk)A{xi)}- (74) 

This provides the counting statistics for an arbitrary multi-terminal struc- 
ture of the type shown in Fig. 4. Below we discuss several examples. 



4.2. NORMAL METAL MULTI-TERMINAL STRUCTURES 

In the case that all terminals are normal metals, we can evaluate Eq. (74) 
further. All terms under the square root are proportional to the unit matrix 
and the trace can be taken easily with the result [45] 

togs 



S(X) = - 



2e 2 



dE x 



(75) 



where gs = J2k=i 9k and t k i = 2g k gi/g^. The argument of the square root 
is the sum over all tunneling events from between the terminals. For two 
terminals we recover the result (44). 

Let us now consider three normal-metal terminals, one of which is op- 
erated as voltage probe, i.e. no average current enters the terminal. This 
layout is depicted in Fig. 4. We assume for the applied potentials that 
V\ < V3 < Vi- Terminal 3 is operated as a voltage probe and it follows, 
that V3 = (gi V\ + 52V2) /(<7i +52)- At zero temperature the energy integra- 
tion yields 

togs V 



s(x) 



+ 



2e 
9i 



92 



1 + t 13 (e i( -X3-xi) - 1) + t^ie^-xi) - l) 



91 +92 



.91 +92 

\J 1 + t 23 (e i ^^ - 1) + ^(etoa-xi) - 1) 



(76) 
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S(x) = (tnie^-^ - 1) + t 23 (e^-^ - 1)) , (78) 



The CGF separates into two terms. The first term corresponds to the en- 
ergy window eVz > E > eVi, in which transport is only possible between 
terminals 1 and 3 or between 1 and 2. The second term results from the en- 
ergy window eV 2 > E > eVJj, in which no electrons can enter into terminal 
2. " 

We now consider a different configuration: an instreaming current is 
divided into two outgoing currents. This corresponds to the voltage config- 
uration V = V\ = V 2 > V3 = 0. The corresponding CGF is 

S(X) = -^Vl + h 3 (e^i-x 3 ) - l) + t 23 ( e ^(x 2 -X3) _ 1) . (77) 

In the limit g 2 , ffi 3> 53 or vice versa the CGF takes the form 

^offs V 
4e 

and the corresponding counting statistics is 

w>) ww- m 

Here iVj = {t^V j e)gigj,j \g\. The statistics is a product of two Poisson 
distributions, i. e. the two transport processes are uncorrelated. 

4.3. SUPERCONDUCTING MULTI-TERMINAL STRUCTURES 

Let us now consider the beam splitter configuration, if the incoming current 
originates from a superconductor. Here, we have to use the NambuxKeldysh 
matrix structure. The layout is as shown in Fig. 4. We choose terminal 3 
as superconducting terminal with V3 = and the potential in the two nor- 
mal terminal is assumed to be the same, V\ = V 2 = V . We consider the 
limit T <C eV <C Exh,A. Transport occurs then only via Andreev reflec- 
tion, since no quasiparticles in the superconductor are present. The Green's 
functions for the terminals can be found in Appendix B.4. We note, that 
the pair breaking effect due to a magnetic field in a chaotic dot with the 
same terminal configuration was studied by Samuelsson and Biittiker [19]. 

We now evaluate Eq. (74) for the three-terminal setup. The CGF de- 
pends only on the differences \i ~ X3 an d X2 — X3i which is a consequence 
of charge conservation and allows to drop the explicit dependence on \3 
below. Introducing pi = 2g 3 gi/ (gj + (gi + g 2 ) 2 ) we find 



S(xi,X2) = 1 x (80) 



\j 1 + \fl+p\ (e i2 » - 1) +p2( e j2 X 2 - 1) + 2pip 2 (e^ 1+X2 ) - ^ 



)• 
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The inner argument contains counting factors for the different possible pro- 
cesses. A term exp(i(xfc + X/) ~~ 1) corresponds to an event in which two 
charges leave the superconducting terminal and one charge is counted in 
terminal k and one charge in terminal /. The prefactors are related to the 
corresponding probabilities. For instance, p\ is proportional to the probabil- 
ity of a coherent tunneling event of an electron from the superconductor into 
terminal 1. A coherent pair-tunneling process is therefore weighted with p\. 
This is accompanied by counting factors which describe either the tunnel- 
ing of two electrons into terminal 1(2) (counting factor exp(«(2xi(2))) — 1)) 
or tunneling into different terminals (counting factor exp(i(xi + \2) — !))• 
The nested square-root functions show that these different processes are 
non-separable. 

It is interesting to consider the limiting case if gz/(g\ + 52) is not close 
to 1. Then, pi^ <C 1 and we can expand Eq. (80) in p\^- The CGF can be 
written as 



The CGF is composed of three different terms, corresponding to a charge 
transfer event of 2e either into terminal 1 or terminal 2 (the first two terms 
in the bracket) or separate charge transfer events into terminals 1 and 2. 
The same form of the CGF appears if the proximity effect is destroyed by 
other means, e. g. a magnetic field [19]. According to the general principles 
of statistics, sums of CGFs of independent statistical processes are additive. 
Therefore, the CGF (81) is a sum of CGFs of independent Poisson processes. 
The total probability distribution P(Aq, A^) corresponding to Eq. (81) can 
be found. It vanishes for odd values of (Ni + N2) and for even values it is 
given by 



Here we have defined the average number of transfered electrons N = 
{toV/e)(g 1 +g 2 ) 2 gl/{{gi+g2) 2 +glf /2 and the probabilities T 1(2) = 51(2)/ 
(91 +52) that one electron leaves the island into terminal 1(2). If one would 
not distinguish electrons in terminals 1 and 2, the charge counting distri- 
bution can be obtained from Eq. (81) by setting \i = X2 = X an d perform- 
ing the integration. This leads to P t s ot (N) = exp(-iV /2)(N /2) N ' 2 /(N/2)\, 
which corresponds to a Poisson distribution of an uncorrelated transfer of 




(81) 




N 1 +N 2 



(82) 
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electron pairs. The full distribution Eq. (82) is given by Pf ot (Ni + N 2 ), 
multiplied with a partitioning factor, which corresponds to the number of 
ways to distribute Ni + N 2 identical electrons among the terminals 1 and 
2, with respective probabilities T\ and T 2 . Note, that T\ + T 2 = 1, since the 
electrons have no other possibility to leave the island. 

5. Concluding Remarks 

Full Counting Statistics is a new fundamental concept in mesoscopic elec- 
tron transport. The knowledge of the full probability distribution of trans- 
fered charges completes the information on the transport mechanisms. In 
fact, the FCS represents all information, which can be gained from charge 
counting in a transport process - a clear progress in our understanding of 
quantum transport. 

In this article we have reviewed the state of the field with particu- 
lar emphasis on superconductor-normal metal heterostructures. We have 
introduced a theoretical method, which combines full counting statistics 
with the powerful Keldysh Green's function technique. This method al- 
lows to obtain the FCS for a large variety of mesoscopic systems. For a 
two-terminal structure a general relation can be derived which contains the 
FCS of all kinds of constrictions between normal metals and superconduc- 
tors. Our method is readily applicable to multi-terminal structures and we 
have discussed one example. An important advantage of the method is that 
it allows a direct numerical implementation, which means that we are able 
to find the FCS of arbitrary mesoscopic SN-structures. 

I would like to thank the Alexander von Humboldt-Stiftung, the Dutch 
FOM, the Swiss SNF, and the NCCR Nanoscience for financial support in 
different stages of this work. Special thanks go to Yuli V. Nazarov, who 
introduced many of the concepts discussed here. Valuable insights emerged 
from discussions with D. Bagrets, J. Borlin, C. Bruder, M. Biittiker, M. 
Kindermann, P. Samuelsson, and F. Taddei. 

A. Field Theoretical Methods 

We summarize some methods and definitions of quantum field theory. In 
the first part we review briefly the standard Keldysh-Green's function tech- 
nique. We follow essentially the review [29]. In the second part we explicitly 
perform the linked cluster expansion for the cumulant generating function, 
which establishes the relation between our definition of FCS and the Green's 
function method. 
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A.l. KELDYSH GREEN'S FUNCTIONS 



One commonly used method to study nonequilibrium phenomena is the 
so-called Keldysh technique. Quite generally time-dependent problems are 
cast into the form of calculating expectation values of some operator A of 
the form (A(t)), where A{t) is the time-dependent operator in the Heisen- 
berg picture with respect to some Hamiltonian: A(t) = U\t,t')A(t')U(t,t'). 
In this expression both the time-ordered evolution operator U(t, t') = T 
exp(— i J*, <ItH(t)) and the anti-time-ordered evolution operator U\t, t') = 

Texp(— i fl cLtH{t)) appear. A diagrammatic theory requires to account 
for the various time orderings, which is rather complicated. A considerable 
simplification arises from Keldysh's trick. We introduce two time coordi- 
nates (ti,i2)j which live on the upper and lower part of the contour (C\, C 2 ), 
and an ordering prescription along the closed time path Ck, depicted in 
Fig. 5. 




Figure 5. The Keldysh time ordering contour Ck 



In the context of full counting statistics we introduce different Hamil- 
tonians #1(2) f° r the two parts of the contour Cu 2 y The actual different 
parts are given by Eq. (13). The rest of the Hamiltonian H sys coincides on 
both parts, just as in the usual formulation [29]. 

We define a contour- ordered Green's function (7~k denotes ordering 
along the Keldysh-contour) 

G CK (t,t') = -i(T K *(t)^(t')). (83) 

Here the field operators can have multiple components, such as spin or 
Nambu for example. This Green's function can be mapped onto a matrix 
space, the so-called Keldysh space, by considering the time coordinates on 
the upper and lower part of the Keldysh contour as formally independent 
variables 

ra +'\ - ( Gn (*'*') Gi 2 (t,t') \ , , 

G(M) = {-G 21 (t,t>) -G 22 (t,t>) )' (84) 

which reads in terms of the field operators in the Heisenberg picture 
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The current is obtained from the Green's functions (with spatial coordinates 
reinserted) as 

j(x,t) = lim (V x - V x 0Gi2(x,t;x',t) (86) 

ZVfl x^x' 

= --i- lim (V x - V x O(Gn(x,t;x / ,t) + G 22 (x,t;x , ,t)), 

Z?7l x^x' 

where the second form is the one we have used in the context of counting 
statistics. The current has still a matrix structure in the subspace of the 
components of the field operators. How the electric current is obtained, 
depends on the definition of the subspace. 

In the usual Keldysh-technique H\ = H2 and we have the general prop- 
erty G11 — G22 = G12 — G21. One element of the Green's function (84) can 
be eliminated by the transformation G = CGC\ where C = (1 — ZT2)/\/2 
and -fj denote Pauli matrices in Keldysh space. Then the matrix Green's 
function takes the form 



G(t,t') = 




where 

G^)(M') = Tie(±(t-mm),&tf)]), (88) 
G K (t,t') = -i({*(t),tf {*)})■ (89) 

The bulk solutions for normal metals and superconductors in App. B.4 are 
given in this form. 

A. 2. THE RELATION BETWEEN CUMULANT GENERATING FUNCTION 
AND GREEN'S FUNCTION 

We establish the connection between the diagrammatic expansions of the 
Green's function, defined in (14), and the cumulant generating function 
in (12). We first show how the diagrammatic expansion of the CGF is 
obtained. In fact, the same expansion occurs in the expression for the ther- 
modynamic potential [46]. According to the definition of the CGF (12) we 
write 

e S(x) = {TcKe -^Ic K d ^ h { q 0) 

We consider a term of the order n in the expansion of the exponent. Such 
a term has a form 

A (~ / • • • / d *l • • • dtn{Tc K h{tl)I c (t 2 ) ■ ■ ■ I c {tn)) ■ (91) 
n\ V 2e/ Jc K Jc K 
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We abbreviate for the sum of connected diagrams of order n 



Qn = ~J ■ ■ ■ [ dh ■ ■ ■ dt n {Tc K Ic(tl)I c {t 2 ) ■ ■ ■ h{tn))c 
n\ Jc K JC K 



(92) 



To count the possible diagrams, which contain p\ connected diagrams of 
order mi, p 2 of order 771,2, and so on, where p\mi + p 2 m 2 + ■ ■ ■ +Pk^k = n , 
we note that their number is equivalent to the number of possibilities to 
assign n operators to p\ cells containing mi places, p 2 cells containing m 2 
places, .... This number is given by 



pi!(mi!)Pip 2 !(m 2 !)P 2 • • -pfc!(m fc !)P* 
and it follows that the CGF can be written as 



(93) 



S(x) = J- — 



P1.P2 

exp 



pi! L 2e 



pi 1 
P2! 



V2 



(94) 
(95) 



Thus, we find that the CGF S(x) is directly given by the sum 

S{x) = Y,{ifj n Qn (96) 

over the connected diagrams only. 

The Green's function introduced in Eq. (14) has a perturbation expan- 
sion 



G(^t; X ',t', X ) = -^EAf-vV/ •••/ dh...dt n x (97) 

^n! \ 2eJ J Ck Jc K 

(r Cjf *(x,t)*t( x ^')/ c ( fl )7 c (t 2 ) . . . I c (t n )) con , 



recalling that all disconnected diagrams are canceled from this expression. 
Now we calculate the current from this Green's function with the same 
current operator used in I c . We find 



i(x,t) 



-„! orhji HTL -L*> x(98) 

(Tc K i c (h)i c {t 2 ) ■ ■ • / c (t„-i)/ c (t)) co „ . 
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In a static situation (as we consider) this current does not depend on t and 
we integrate along the Keldysh contour between and to- Using Eq. (92) 
it follows that 

oo , y \n-l 

2t /(x) = -E n (-^j Qn (99) 

n=l 

By comparing the right-hand sides of (99) and (96) we finally obtain the 
relation between the x-dependent current and the CGF: 

7 (x )=if^. (100) 
to ox 

The constant contribution to S(x) follows from the normalization 5(0) = 0. 



B. Quasiclassical Approximation 

In practice an exact calculation of Green's functions is impossible in virtu- 
ally all mesoscopic transport problems. An important simplification is the 
quasiclassical approximation [30] , which makes use of the smallness of the 
most energy scales involved in transport with respect to the Fermi energy 
Ep. We briefly summarize the derivation of the basic equations [29]. Here 
we concentrate on the derivation in the context of superconductivity. The 
inclusion of spin-dependent phenomena is straightforward. 



B.l. EILENBERGER EQUATION 

The starting point is the equation of motion for the real-time single-particle 
Green's function. We consider here the static case, in which the equations 
can be considered in the energy representation. The equation of motion 
reads 

~2 



Ed?. - 



2m 



+ Ep — £ 



G(x,x';£7) =<5(x-x') 



(101) 



Here the " denotes matrices in the combined Nambu x Keldysh-space and we 
use <7j(fj) to denote Pauli matrices in Nambu(Keldysh)-space. Ep = p F /2m 
is the Fermi energy and the selfenergy X includes scattering processes. The 
complicated dependence on two spatial coordinates in Eq. (101) can be 
eliminated by the following procedure. We introduce the Wigner transform 

G(r,p;E)= J d 3 s exp (ips) G(r + s/2, r - s/2; E) , (102) 

for which the equation of motion reads 



Ea 3 - ivV r - E(r, p) G(r, p,E) = l. 



(103) 



31 



In this equation, we can neglect the dependence on the absolute value of the 
momentum in the expression in the brackets, since G(r, p, E) is strongly- 
peaked at p = pp. We now subtract the inverse equation and integrate the 
resulting equation over £ = p 2 /2m. We are lead to the definition of the 
quasiclassical Green's function 

g(r,p F ,E) = - I *™G(r,p;E). (104) 
7r J m 

The quasiclassical Green's function obeys the Eilenberger equation [30] 

^-Vj(r,v F ,£) = [-iEa 3 + ia{r,v F ),g{r,v F ;E)] . (105) 

The current density is obtained from the generalized matrix-current 

j(r,v F ,£) = e 2 iVov F 5(r,v F ,£) (106) 
j(r) = j-J dETrf K (](r,v F ,E)) PF , (107) 

where () PF denotes angular averaging of the momentum direction and fx = 

Physically, we have in this way integrated out the fast spatial oscilla- 
tion on a scale of the Fermi wave length A F , and the remaining function g 
is slowly varying on this scale. We have also assumed that the self energy 
S does not depend on the momentum p. This has two important conse- 
quences. First, the equation is now much simpler and has the intuitive 
form of a transport equation along classical trajectories. It is in fact easy 
to show, that Eq. (105) reproduces the well known Boltzmann equation 
for normal transport. Second, we have a price to pay. Obviously the above 
derivation fails in the vicinity of interfaces or singular points. This means, 
that the boundary conditions have to be derived using the underlying non- 
quasiclassical theory. Fortunately, this can be done quite generally and one 
simply has to use effective boundary conditions. Related to the problem of 
boundary conditions is the homogeneity of Eq. (105). In fact, it was shown 
by Shelankov [47] that the condition ensuring the regularity for x — ► x' in 
the general Green's function leads to a normalization condition 

5 2 (r,v F ,£) = l. (108) 

This condition is very important, since almost any further manipulation of 
Eq. (105) relies on it. 

In the context of superconductivity the most important contribution to 
the self energy is the pairpotential 

& R(A) = _,A , A = -i\ J dE(g« di J PF , (109) 
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where 'offdiag' denotes that only the off-diagonal components in Nambu 
space should be considered. A is the attractive BCS interaction constant. 

B.2. THE DIRTY LIMIT - USADEL EQUATION 

An important simplification arises if the system is almost homogeneous in 
the momentum direction. This is e.g. the case for diffusive systems, in which 
the dominant term in the self-energy arises from impurity scattering. In 
(self-consistent) Born approximation the impurity self-energy has the form 

a(r,E)=--^(g(r,p F ,E)) PF . (110) 

^ ' imp 

In the limit l/rj mp S> E,A, etc., the Green's functions will be nearly 
isotropic and we make the Ansatz g(r, pp, E) = G(r, ^)+ppG(r, E). Using 
Eq. (105) together with the normalization condition we obtain the so-called 
Usadel equation [31] 

Vj(r,£) = \-iE + A,G(r,E)] , (111) 



e 2 N 

where the generalized matrix current is now 

'j(r,E) = aG(r,E)VG(r,E) (112) 

The conductivity is given by the Einstein relation a = e 2 NoD, with the 
diffusion coefficient D = VpTi mp /3. The current density is 

3(r) = ±-JdETr[f K ](r,E)]. (113) 



B.3. BOUNDARY CONDITIONS 

Close to boundaries the quasiclassical equations are invalid and have to be 
supplemented by boundary conditions. In the general case these boundary 
conditions have been derived by Zaitsev [48]. They are rather complicated 
and we will not treat them here. In diffusive systems a very concise form 
of the boundary conditions was obtained by Nazarov[33]. Under the as- 
sumption that two diffusive pieces of metals are connected by a quantum 
scatterer, the matrix current depends only on the ensemble of transmission 
eigenvalues and has the form (21). The boundary condition is equivalent to 
the conservation of matrix currents in the adjacent metals. 
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B.4. BULK SOLUTIONS 



We summarize the necessary ingredients for the various circuits treated in 
this review. As terminals we consider only normal metals or superconduc- 
tors. They are determined by external parameters like applied potentials or 
temperature. The matrix structures are obtained from the bulk solutions 
of the Eilenberger- or Usadel-equations. We give below their form in the 
triangular Keldysh- matrix representation (87). 

A normal metal at chemical potential fj, and temperature T e is described 
by a Green's function 

G N (E) = a 3 f 3 +(f 1 +if 2 ) (1 - f(E) - f(-E) + a 3 (f(-E) - f(E))) (114) 

with the Fermi distribution f(E) = (exp((£ - fi)/k B T e ) + 

A superconducting terminal at chemical potential (is = is described 

by 



G S (E) = 1 (R + i) + X - (R - i) (f 3 + (fi + if 2 ) tanh ( 



E 



2 V J 2 V 7 V \2k B T e 

(115) 

where the retarded and advanced functions are 

r(A) = (7 a) fR(A) ) = AE±m+2L, (116) 

V JR{A) ~9r(A) J ^E±i5) 2 - |A| 2 

where 5 is a broadening parameter. The gap matrix contains the dependence 
on the superconducting phase (f> 

X ( |A|e^ \ 
= ( |A|e^ J • (U7) 

In the limit of zero temperature and \E\ <C A the bulk solutions simplify 
to 

g n (e) = {^-f}+ if }\ s ^ eV ( i, ' !f'<! e !;! (us) 

G s ((f)) = u\ cos(<^>) — o"2 sin(<^>) . (119) 



C. Appendix: CGF for the Single Channel Contact 

C.l. SUPER-NORMAL CONTACT: COEFFICIENTS A N 

We present the coefficients in the CGF (37) for a channel of transparency 
T n . We assume the superconductor to be in equilibrium and the normal 
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metal at a chemical potential /j,n. The occupation factors are f± = fo(±E— 
Hn) and /± = fo(±E), where fo is the Fermi-Dirac distribution. The coef- 
ficients take the form 



K,i = T n (l-T n /2) 



A 



n,2 



2(#R - 9 a) 



+2T, 



2-T n (g R -l)){2-T n (g A + l)) 
/£(l-/£)+/5(l-/£) 
1 - JrJa - 9R.9A 



(120) 



n (2-T n {g R -l))(2-T{g A + l)) 

(tf - f N )U s + - f s )(i - (tf - /^)(/| - Z 5 )) 

+2(/|-/f) 2 (l-/^)(l-/i v )" , 



T^2 

2 



(121) 



1 + Ma ~ 9R9A ~ (f+ ~ /_r(l - 7fl/A - 9R9A) 
(2-T n (g R -l))(2-T n (g A + l)) 



The coefficients A n _ q can be obtained from Eq. (120) and Eq. (121) by 
the substitution /+ ^ <-> 



(1 — f^ N ^), i.e. interchanging electron-like and 



hole-like quasiparticles. 



C.2. SUPER-SUPER CONTACT: COEFFICIENTS A 



N 



Introducing q = (1 — <7r<7a)(1 ~ ^- 2 ) + /r/a(1 + ^- 2 ) the coefficients may be 
written as 



A n,±2 



Qn 



16 



T 2 

64 y ' 
T n 



l-r n /^sin^(| 



<? - 4/r/a sin 2 



(/r + /a)/icos - T - /a) sin - 



1- T^/isin 2 (|) 



(122) 
(123) 

(124) 
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